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THE  MOUNTAIN  IRON 
DIFFUSION  PROGRAM:  PHASE  I 

W  T.  Hinds  and  P„  W.  Nickola 

ABSTRACT 

Field  diffusion  experiments  were  conducted  at  Vandenberg  Air 
Force  Base,  California,  during  1965  and  1966.  This  program,  nicknamed 
"Mountain  Iron,  "  was  undertaken  to  establish  quantitative  diffusion  pre* 
dictions  for  use  as  range  safety  tools  in  the  "South  Vandenberg"  ballistic 
and  space  vehicle  operations,  The  only  pollutant  source  character  studied 
was  toxic  propellant,  materials  released  continuously  from  ground  level 
pools. 

Volume  I  of  this  report  contains  the  operational  applications  and 
limitations  of  the  equations  and  other  results  obtained. ^Volume  II 
includes  a  description  of  the  conduct  of  diffusion  experiments,  discussion 
of  results,  summaries  of  experimental  data,  and  an  analysis  of  the 
collected  data,  including  development  of  the  various  diffusion  equations, 
preliminary  climatology,  and  a  discussion  of  applicability  limitations  and 
recommendations  for  further  studies.,  A  third  item  of  the  report  is  a 
compilation  of  experimental  data  obtained.  These  data,  generally 
described  m  Volume  II,  are  too  voluminous  to  be  distributed  to  all  recipi  ¬ 
ents  of  Volumes  I  and  II.  However,  the  data  are  archived  at  the  Environ¬ 
mental  Technical  Applications  Center  (ETAC),  1210+h  Wea  Sq,  Bldg.  159 
Navy  Yard  Annex,  Washington,  D.C.  20370,  from  which  authorized 
agencies  may  obtain  the  data  through  accepted  channels  of  request. 
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THE  MOUNTAIN  IRON 
DIFFUSION  PROGRAM:  PHASE  I 

PREFACE 


The  prime  considerations  of  Volume  I  are: 

•  Compliance  with  the  Clean  Air  Act,  Public  Law  88-206,  December, 
1963 

•  Compliance  with  AFR  161-22,  Environmental  Pollution  control, 
as  supplemented 

•  Protection  of  life  and  property  through  adequate  range  safety 
procedures. 

The  details  and  basic  data  pertinent  to  the  final  equation  formulation  will 
be  found  in  Volume  II;  the  descriptions,  statements  of  applicability,  usage 
and  limitations  of  results  are  discussed  in  this  volume. 
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THE  MOUNTAIN  IRON 
DIFFUSION  PROGRAM:  PHASE  I 

W.  T.  Hinds  and  P.  W.  Nickola 

INTRODUCTION 

The  objective  of  Project  Mountain  Iron  -  Phase  I  South  Vandenberg 
was  to  develop  a  system  for  describing  the  potential  of  the  atmosphere  for 
diffusing  and  diluting  pollutants  employed  as  propellants  in  launch  prepar¬ 
ation  and  operations  for  the  South  Vandenberg  area. 

Mesoscale  atmospheric  diffusion  parameters  are  introduced  into  range 
safety  considerations  through  a  computer  controlled  data  acquisition  and 
display  system.  Diffusion  equations  are  either  applied  by  computer  pro¬ 
gram  to  predict  final  extent  of  the  existing  potential  toxic  hazard  area,  or 
by  forecaster  manipulation  through  tables  or  computer  program. 

Preliminary  thought  on  the  subject  indicated  that  diffusion  equations 
valid  for  one  area  of  the  range  would  not  be  valid  for  another  (as  this 
chapter  will  indicate)  because  of  terrain  and  weather  regime  differences. 

Present  instrumentation  for  the  collection  of  measurable  parameters 
was  proven  to  be  adequate  for  insertion  into  the  diffusion  equations  for  pre¬ 
diction  of  potential  or  actual  toxic  hazard  areas  at  the  95%  confidence  level 
as  a  minimum. 
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CHAPTER  1. 

SOUTH  VANDENBERG  EQUATIONS 

The  equations*  presented  here  are  in  the  forms  utilized  in  tables 
and  computer  programs  related  to  the  AFWTR  meso-meteorologioal 
data  processing  system  with  application  to  prediction  of  atmospheric 
difl’usivity,  A  discussion  of  the  limitations  on  the  use  of  the  equations 
is  included  to  assure  judicious  use  of  the  equations  and  predictions. 

DERIVATION  OF  THE  EQUATIONS 

Two  parallel  approaches  were  pursued  during  the  analysis  of  the 
Mountain  Iron  data:  one  was  base'd  upon  physical  reasoning  and  using  time 
as  the  downwind  variable;  the  other  was  a  purely  statistical  approach  with 
distance  as  the  downwind  variable.  The  difference  between  the  two  equa¬ 
tions  is  minimal,  but  the  statistical  approach  yielded  slightly  more 
accurate  predictions  and  substantially  si-  ••  er  data  manipulations. 

RECOMMENDED  AND  ALTERNATE  EQUATIONS 

In  the  development  of  the  statistical  model,  two  equations  of  almost 
identical  format  were  finally  chosen  to  provide  the  best  available  pre¬ 
dictions.  The  recommended  equation  uses  the  temperature  difference 
between  C>  and  200  feet;  the  deeper  layer  parameter  produces  a  slightly 
improved  reliability  of  prediction  within  a  factor  of  2  when  compared  to  the 
results  from  a  fi  to  54  foot  temperature  difference. 

At  the  95%  level,  the  distance  to  a  given  concentration  is 

X  =  790(X/Q)’0,  °4  Cg0,  192  if  470  (AT200  +  10.  S)°*  6  H)  (1) 

This  equation  refers  to  average  concentration  of  NO,^  only,  in  parts  per 
million  The  release  rate,  which  was  used  to  normalize  the  concentration, 
is  units  of  pounds  per  minute. 


A  ilet ailed  discussion  of  the  derivation  of  prediction  equations  is  given  in 
Chapter  (i,  Volume  II. 
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Because  6  to  300  foot  temperature  differences  are  available  at  very 
few  sites  on  South  Vandenberg,  an  alternate  equation,  developed  for  6  to 
54  foot  temperature  difference,  is  used  here  to  form  an  equation  for 
operational  use  in  the  same  manner  as  the  recommended  equation.  At 
the  95%  level,  the  distance  equation  is 


X  =  608(X)  -  55 


-0.229  n-0.566  +  9)0.  852 

b  6 


(2) 


Although  this  alternate  equation  differs  slightly  from  the  recommended 
equation,  it  still  refers  to  parts  per  million  of  NC>2  only,  and  is  normalized 
to  a  release  rate  in  terms  of  pounds  per  minute. 


EQUATIONS  FOR  AN  ARBITRARY  CONTAMINANT 

If  trouble  from  some  other  gas  is  expected,  the  equations  must 
be  corrected  to  allow  for  variation  in  molecular  weight.  If  M  is  the 
molecular  weight  of  the  gas  of  concern,  then  the  95%  confidence  level 
equation  in  terms  of  distance  to  a  concentration  is 


X  =  5000(^) 


-0.55 


o-0.229  u-0.566  (AT54  + 


9) 


0.852 


(3) 


If  the  contaminant  is  a  particulate  rather  than  a  gas,  specification 
of  parts  per  million  loses  its  significance;  that  is,  the  mass  of  the  con¬ 
taminant  per  cubic  meter  is  the  pertinent  measure.  At  the  95%  confidence 
level,  with  X  in  grams  per  cubic  meter  and  Q  in  pounds  per  minute,  the 
distance  to  a  concentration  is 


X  =  19.  8(h~°‘  55  op0,  229  u'°*  566  ^T^4  +  9)°’  852 


(4) 


TABLES  AND  COMPUTER  PROGRAMS  FOR  TOXIC  HAZARD  AREA 
COMPUTATION 

Since  there  is  no  provision  for  machine  display  of  a  forecast  planform 
of  a  potential  or  actual  toxic  hazard  area,  some  means  of  rapid  computation 
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must  be  made  available  to  the  responsible  agency.  Equation  (6)  page  8, 
Chapter  2,  forms  the  basis  for  computer  prepared  tables  for  South 
Vandenberg.  These  tables  enable  the  finding  of  X  when  any  valueB  of  u , 
A,  Ogand  AT  are  given. 

The  computer  program,  utilized  in  the  mesoscale,  real  time, 
meteorological  parameter  sampling  system,  presents  displays  of  u,  Og, 

AT  and  wind  direction  together  with  X  for  various  A  values  for  current, 
potential,  toxic  hazard  area  description 
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The  equations  presented  in  Chapter  1  have  been  manipulated  into 
other  forms  from  which  tables*  and  computer  programs  have  been 
prepared. 

All  equations  in  this  chapter  are  based  upon  a  ground  pool  of 
(or  NOg)  that  provides  a  continuous  point  source  of  pollutant,  The 
equations  will  provide  at  least  95%  confidence  that  the  actual  or  potential 
toxic  hazard  areas  defined  by  the  equations  will  contain  all  concentrations 
of  pollutants  equal  to  or  in  excess  of  25  ppm. 

The  South  Vandenberg  Diffusion  Equation  is  a  manipulation  of 
Equation  (2)  which  becomes 


X  «  103  P  a 


-0.229  /a~54 


(ATg  +  9) 


0.  852 


where  X  is  distance  in  feet  from  source  to  where  the  concentration  of 
pollutant  becomes  less  than  25  ppm, 

o„  is  the  standard  deviation  of  wind  direction  in  degrees  at  the 


ATg  is  the  temperature  difference  in  degrees  Fahrenheit  between 
the  heights  of  6  and  54  feet  above  ground  level. 

P  is  a  modified  release  rate,  calculated  from 


P  *  0.  162  u'0,  126  A0,  550 


where  u  is  the  average  wind  speed  in  knots  at  12  feet  above  ground  level 


Tables  for  finding  P  =  f  (u,  A)  and  X  =  f(P,  Oq,  T)  were  prepared  for  use 
on  South  Vandenberg  AFB. 
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and  A  is  the  area  in  square  feet  of  the  ground  level  pool  of  pollutant.  Thus, 
combining  Liquations  (8)  and  (5a), 


,  r.n/.  _-0.  128  .  0.  550  -0.229 


X  =  16. 686  u 


(ATg  +  9) 


0.  852 


The  North  Vandenberg  Diffusion  Equation  is  the  result  of  the  Dry 
Gulch  portion  of  The  Ocean  Breeze  and  Dry  Gulch  Diffusion  Programs/1^ 
This  equation  from  which  tables*  and  computer  programs  were  prepared  is 


X  =  9.  82(1.71)( 


^)"°'5°V-25e  <4T|4  +  t0>a-21 


where  Q  is  the  release  rate  in  pounds  per  minute  calculated  from 
Q  =  0. 0242  u0, 8  A 

with  IT  the  average  wind  speed  in  feet  per  second  at  12  feet  above  ground 
level  and  other  terms  are  as  for  Equation  (6).  Through  algebraic  manipu¬ 
lation  and  conversion  of  u  to  knots, 


X  =  0. 601  u 


_0. 4072  .  0.509  -0.258,.  -.54 


(ATjj  +  10)‘ 


A  comparison  of  Equations  (6)  and  (8)  reveals  that  similar 
measurable  meteorological  parameters  are  involved  in  the  South  Vandenberg 
and  North  Vandenberg  diffusion  equations.  Less  apparent  are  the  effects 
of  these  parameters  on  X.  Since  the  equations  are,  in  essence,  empirically 
derived  from  prior  studies  and  theory  by  statistical-fit,  some  variance  is 
found  in  the  value  assigned  to  coefficients  and  exponents  when  the  North  and 
South  Vandenberg  equations  are  compared. 


*  Tables  for  finding  Q  =  f(u,  A)  and  X  =  f(Q,  o,,.  AT)  were  prepared  in  a 
previous  study  for  use  on  North  Vandenberg  AFB. 
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The  North  Vandenberg  equation  shows  a  systematic  error  when 
applied  to  South  Vandenberg  data  (Figure  1).  On  the  average,  the 
equation  overpredicts  by  about  a  factor  of  2  but  underpredicts  26%  of  the 
cases  examined.  The  relative  scatter  increases  as  exposures  decrease 
(at  greater  distances).  Because  of  limitations  in  sampling  grids,  com¬ 
parisons  are  taken  out  only  to  8  kilometers. 

Generally,  the  predictions  from  the  North  Vandenberg  equation 
are  rather  good.  The  major  drawback  to  its  use  is  the  too  rapid  decrease 
of  exposure  with  distance,  which  would  lead  to  underprediction  if  pressed 
beyond  several  kilometers.  Wind  speed  was  found  not  to  be  an  influ¬ 
ential  parameter  in  the  North  Vandenberg  equation,  which  is  also  too 
sensitive  to  temperature  stratification  to  be  accurate  over  South 
Vandenberg. 

Sensitivity  of  Equations  to  Parameters 

Comparing  Equations  (6)  and  (8),  term  by  term,  fob'll: 

North  X  =  f(u) 

South  X  =  f(l/u) 

This  reflects  the  greater  turbulent  flow  and  mixing  over  South  Vandenberg 
due  to  the  rougher  terrain.  Given  the  same  u  -  10  kts: 

North  u0, 4072  =  2.  554 
South  u'0, 126  «  0.7482 

The  u  effect  for  South  Vandenberg  is  much  more  pronounced  (at  least  by 
a  factor  of  three)  than  for  North  Vandenberg.  For  A: 

North  X  =  f(A) 

X  =  f(A> 


South 
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Predicted,  E^/Q^. 


FIGURE  1.  Mountain  Iron  Data  Versus  Predictions  from  Weather 
Information  Network  and  Display  (WIND) 


A.  A 


r$i\wi,-:i72  voj.  ] 


Given  the  same  A  -  10,000  square  feet: 


North 

South 


A°'®09  =  108.65 
A0'550  *  105.20 


For  AT®4: 

D 


North 

South 


X  =  f(ATg4) 
X  =  f(AT®4) 


Given  the  same  AT®4  =  1. 0  °F: 


North  (AT®4  +  10)2,21  =  169.0 
South  (AT®4  +  9)0,852  s  7<  j 


For  a 


6: 


North 

South 


x  ‘  «£> 


X  - 


f(~) 

le 


Given  the  same  o^  =  13  “F; 


North 

South 


-0.258  =  1.938 

,"0. 229  =  1. 799 

'9 


Substituting  these  values  in  Equations  (1)  and  (8): 

U  A  Cy 

North 


T 


South 


X=  0.601  (2.554)  (108.65)  (1.938)  ( 1 69.  9)  -  54 ,  900  f0t., 
X=  16.86  (0.7582)  (105.20)  (1.799)  (  7.  1)  -  1 7,  200  feet 
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r>  1  n  VINO  V.  VIN 


nt-J  USB  OK  THB  SOUTH  VANDENEERG  EQUATION (S) 


Results  from  any  empirical  study  arc  limited  in  their  usefulness, 
and  the  statistical-fit  equations  developed  from  the  Mountain  Iron  data  are 
ih>  exceptions.  To  offer  some  guidance  in  the  judicious  use  of  the  equa¬ 
tions,  the  limits  imposed  by  test  conditions  are  briefly  examined  here. 


Source*  I .ocution  Limitations 


The  equation  specifically  applies  to  the  VIP-1  area.  However,  it 
should  he  applicable  to  source  points  resembling  the  Source  A  area:  near 
the  ocean,  with  unobstructed  upwind  fetch,  and  with  a  fairly  rough  downwind 
region  over  which  diffusion  takes  place  (which  assures  a  level  of  mixing 
more  or  less  comparable  to  that  observed  from  Source  A).  These  restric¬ 
tions  mean  that  the  equation  must  not  be  expected  to  apply  to  diffusion  from 
a  source  in  Honda  Canyon ,  for  example,  or  in  the  Santa  Ynoz  Valley ,  especially 
in  westerly  winds.  Likewise,  a  source  in  the  Sudden  Ranch  area  can  not 
he  reliably  simulated  from  the  Phase  I  Mountain  Iron  equation.  Source 
points  that  are  reasonable  extrapolations  from  the  Source  A  area  lie 
he  tween  Pale  I,  Pale  II,  Pad  A  and  Pad  D--a  fairly  generous  portion  of 
till'  South  Vandcnberg.  Probably,  any  departure  from  the  immediate 
area  of  Source  A  will  result  in  degrading  t ho  accuracy  and  reliability 
of  the  predict  ions.  However,  in  this  respect,  comparisons  of  the 
Source  H  tests  with  predictions  From  the  Source  A  equation  and  tower  300 
data  are  encouraging.  The  Source  B  tests  show  that  exposures  from  this 
release  point  can  be  expected  to  be  comparable  to  those  from  Source  A 
over  moderate  distances. 


Source  Typo  Restrictions 

Most  of  tlu*  Mountain  Iron  tests  were  of  15  or  30  minute  duration ,  with 
a  constant  and  continuous  release  rate,  while  9  tests  were  only  five  minutes. 
Since  no  discernible  difference  exists  between  15  and  30  minute  tests,  these 
were  treated  in  the  analysis  as  a  homogeneous  group.  The  5  minute  tests 
presented  a  slightly  different  result:  they  hiiucated  a  more  rapid  decrease 
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in  exposure  vvit.li  distance  than  the  other  tests  in  agreement  with  theoret¬ 
ical  prcdic  ons  and,  1'or  that  reason,  are  not  strictly  within  the  domain 
of  the  equations  recommended  for  use. 

The  most  severe  restriction  on  source  type  is  the  height  of  release, 
with  only  ground  level  releases  being  pertinent.  If  any  substantial  eleva¬ 
tion  from  ground  level  is  postulated,  t lie  Mountain  Iron  results  are  not 
applicable.  Thus,  bouyant  plumes- -as  from  a  conflagration- -are  among 
the  excluded  types  of  sources;  so  are  simple  elevated  sources.  Some  use 
of  the  equation  can  be  made  as  a  limit,  which  slightly  elevated  sources  will 
soon  approach,  but  there  is  little  possibility  of  accurately  predicting  the 
position  or  magnitude,  of  the  maximum  exposure  from  any  elevated  source. 

It  seems  probable  that  the  equation  will  not  servo  even  as  a  limiting  ease 
if  the  source  is  high  (such  as  a  couple  hundred  feet),  or  if  the  elevated 
release  is  into  a  stable  atmosphere.  Only  a  substantially  modified  testing 
program  will  yield  results  pertaining  to  bouyant  or  elevated  sources.  As 
in  the  case  with  the  North  Vandenberg  WIND  equation,  the  South  Vandenberg 
equation  applies  most  accurately  to  a  cold  spill  situation,  especially  one 
with  a  limited  horizontal  extent  since  the  tests  were  conducted  with 
essentially  point  sources. 

Puffs,  even  if  not  bouyant,  are  not  correctly  handled  by  the 
recommended  equation  because  of  the  great  difference  in  the  mechanism 
governing  diffusion  of  puffs  and  plumes  (the  5-minute  tests  approach 
puffs  in  behavior).  In  particular,  puffs  during  a  portion  of  their  history 
diffuse  faster  than  plumes,  leading  to  somewhat  lower  exposures  for  a 
given  release  rate.  Without  data,  no  guidance  is  possible  other  than 
theoretical  developments  beyond  the  scope  of  this  report. 

Distance  I, innits 

All  Mountain  Iron  tests  were  c.-  fined  to  distances  less  than 
about  10  miles,  with  only  a  very  few  tests  successful  in  providing  data  at 
that  distance.  The  majority  of  the  tests  terminated  at  5  miles  or  so, 
since  the  sampling  was  restricted  to  the  South  Vandenberg  area. 
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For  this  rc.*.-u«i,  confident  use  of  the  equations  must  be  restricted  to  per¬ 
haps  10  or  11  miles.  The  systematic  error  pointed  out  in  the  statistical 
fit  portion  of  Chapter  (1,  Volume  II  means  that  extreme  caution  is  required 
even  at  that  distance  because  underprediction  is  almost  assured  during 
daytime  conditions  beyond  some  8  or  f)  miles.  There  is  no  help  for  this 
unless  a  quite  artificial  correction  is  applied  to  daytime  tests  to  reduce 
the  error,  and  such  a  practice  is  not  recommended.  More  work  would 
he  required,  and  substantially  different  testing  programs  carried  out, 
to  provide  a  real  basis  for  prediction  of  exposures  at  distances  exceeding 
10  miles. 

DISCUSSION  OF  PLUME  GEOGRAPHY 

Plume  Width  with  Regard  to  Atmospheric  Stability 

An  analysis  of  Sutton's^  diffusion  equations  for  a  continuous 
point  source  indicates  that  the  area  enclosed  by  an  isopleth  drawn  through 
concentrations  of  equal  value  will  be  "cigar  shaped,  "  with  the  major  axis 
along  the  direction  of  the  moan  wind  and  the  maximum  cross-wind  dimen¬ 
sion  of  the  isopleth  occuring  about  half  way  between  the  source  and  the  end 
of  the  plume.  Using  Sutton's  theoretical  works,  it  has  been  demonstrated 
and  confirmed  by  experimental  data  that  over  level  terrain,  the  width  of 
the  plume  for  unstable  conditions  is  one-third  the  length;  for  neutral  con¬ 
ditions  one-fifth  the  length;  and  one-eighth  the  length  for  stable  conditions. 
The  angular  width  of  the  plume  is  a.  direct  function  of  time  and  wind  direc¬ 
tion  variability.  ^  Modification  of  these  results  by  flow  over  mountainous 
regions  is  easily  seen  in  the  plume  patterns  of  Appendix  A. 

Effect  of  Terrain 

In  general,  the  plumes  observed  over  South  Vandenberg  can  be 
said  to  fall  into  one  of  three  types; 

•  A  modified  cigar 

•  A  canyon  low 

•  A  multiple  maximum  type. 
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The  modified  cigar  type  is  only  slightly  different  than  the  ordinary  plume 
from  sources  over  flat  terrain;  they  may  be  characterized  by  bunding  or 
turning,  but  represent  a  monotonic  decline  of  exposure  with  distance. 

The  canyon  low  type  results  from  flow  over  two  ridges  with  a  deep,  steep 
canyon  between.  The  exposures  found  in  the  canyon  are  significantly 
lower  than  on  either  ridge.  The  multiple  maximum  type  is  really  an 
elaborate  canyon  low  type,  resulting  from  flow  over  two,  three,  or  even 
more  ridges  and  valleys.  As  might  be  expected  from  a  look  at  the 
terrain  of  South  Vandenberg,  the  field  samples  were  strongly  biased 
toward  the  canyon  low  type,  and  probably  the  canyon  low  type  can  be 
expected  to  occur  most  of  the  time  with  northwest  flow  (the  prevailing 
wind  direction)  from  a  source  near  the  north  side  of  South  Vandenberg. 
Figure  2  shows,  dramatically,  the  clear-cut  character  observed  during 
an  actual  test. 

OPERATIONAL  USE  OF  AIRCRAFT  SAMPLING 

The  Battelle-Northwest  tracer  sampling  aircraft  was  operated 
successfully  during  23  of  the  field  tests  of  the  Mountain  Iron  Diffusion 
Program.  The  guides  and  techniques  suggested  are  from  data  collected 
by  the  aircraft,  and  detailed  discussion  of  aircraft  sampling  is  included 
in  Appendix  B,  Volume  II. 

The  reader  is  reminded  that,  although  the  aircraft  sampling  was 
obviously  done  aloft,  the  tracer  source  in  all  field  tests  was  within  a 
few  meters  of  the  surface. 

Vertical  Extent  of  a  Tracer  Released  Near  the  Surface 

Measurements  of  concentration  made  by  an  aircraft  of  a  tracer 
released  from  a  point  near  the  surface  revealed  the  following: 

•  Where  only  superadiabatic  and/or  adiabatic  conditions  existed, 
the  tracer  was  found  at  relatively  high  elevations  and  at  rela¬ 
tively  low  concentrations.  Under  these  conditions,  tracer  was  found 
as  high  as  900  feet  above  the  surface  at  distances  less  than  2 
kilometers  from  the  source;  2000  l'eet  at  4  kilometers;  and  2000  feet 
at  18  kilometers. 
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FIGURE  2.  Canyon  Low  Type  of 
(MI  Test,  Number  48) 
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•  Where  thermal  stratifications  between  isothermal  and  adiabatic 
existed,  the  tracer  was  frequently  found  well  within  or  at  the  top 
of  these  relatively  stable  layers,  but  heights  were  not  as  great  as 
mentioned  above. 

•  Where  a  temperature  inversion  existed,  tracer  penetration  into 
the  inversion  layer  was  minimal. 

Horizontal  Extent  of  Aloft  Tracer  with  Respect  to  Ground  Plume 

The  measured  or  predicted  ground  level  plume  centerline  can  be 
used  as  a  basis  for  specifying  the  expected  location  of  the  aloft  portion  of 
the  plume.  The  first  step  in  this  procedure  is  to  consider  the  aloft  plume 
at  a  chosen  distance  (irrespective  of  altitude)  as  being  projected  to  the  sur¬ 
face.  Figure  3  then  specifies  the  distance  from  the  ground  plume  center- 
line  where  the  projection  will  fall  50  and  83%  of  the  time.  For 
instance,  the  ground  projection  of  the  aloft  plume  at  8  kilometers  can  be 
expected  to  be  found  within  about  1.2  kilometers  of  the  ground  plume 
centerline  50%  of  the  time,  and  within  2.4  kilometers,  83%  of  the 
time. 

The  relationship  is  based  on  field  measurements  to  distances  of 
10  kilometers  from  the  source.  Beyond  that,  the  plume  centerline  was 
not  defined.  Thus,  the  dotted  portion  of  the  curves  beyond  10  kilometers 
from  the  source  are  extrapolations. 

An  alternate  method  of  estimating  the  location  of  the  aloft  plume 
is  available  if  the  lateral  extent  of  the  ground  plume  is  known  or  is 
predictable.  Approximately  90%  of  the  intercepts  of  the  aircraft 
with  the  tracer  were  demonstrated  to  have  occurred  directly  above  the 
area  swept  by  the  ground  level  plume. 

Comparison  of  Peak  Aloft  Concentration  to  Ground  Concentrations 

A  parameter  of  interest  is  the  highest  concentration  that  can  bo 
expected  in  the  zone  above  a  ground  level  plume.  If  one  knows  or  can 
reasonably  forecast  ground  level  plume  centerline  mean  concentrations, 
then  an  estimate  of  the  peak  instantaneous  concentration  can  be  made 
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for  the  alert  plume.  The  relationship  is  based  on  exposure  data  from 
filters  exposed  at  the  1.5  meter  elevation,  and  on  concentration  measure¬ 
ments  made  at  a  variety  of  elevations  by  the  aircraft. 

A  good  estimate  for  this  relationship  is 

W2-0*1'2  (9) 


where  X.  is  the  near  instantaneous  aloft  peak  concentration  at  distance 
ip 

X,  Tp  is  the  mean  ground  level  centerline  concentration  at  X,  and  X  is  the 
distance  from  the  source  in  kilometers. 


Figure  4  plots  the  ratio  y  /Y  for  distances  out  to  about  26 

ip  P 

kiloni  aters--the  maximum  distance  to  which  data  were  obtained.  The 
dotted  lines  suggest,  the  envelope  ofx.  /"X  values  observed  from  all  field 
tests  where  such  computations  were  possible.  Their  magnitudes  are 
approximately  1/5  and  5  times  the  "typical"  values  shown  by  the  solid 
curve. 


No  instantaneous  measurements  of  tracer  concentration  were 
made  at  the  ground  level,  so  ratios  of  instantaneous  ground  level  peak 
concentrations  to  X,  cannot  be  directly  investigated.  However,  beyond 
distances  of  about  5  kilometers,  it  is  unlikely  that  the  ratio  of  ground 
level  instantaneous  peak-to-mean  centerline  concentration  will  exceed 
the  ratio  whoso  equation  is  given  above. 

SECONDARY  INVESTIGATIONS 

In  conjunction  with  the  work  leading  to  the  prediction  equation- - 
the  object  of  this  study- -several  subsidiary  facets  of  diffusion  over 
South  Vandenberg  were  studied.  Although  generally  brief  and  not  at  all 
definitive,  these  studies  were  certainly  of  interest  in  the  conclusions 
that  were  possible.  The  topics  included  are: 

•  Short-term  releases 

•  Trajectory  analysis  of  some  of  the  tests 

•  Ways  of  constructing  trajectories. 
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The  initial  effort  made  here  on  trajectory  studies  indicates  that  investi¬ 
gation  of  transport  on  a  rnesoscale  could  be  particularly  promising  for 
future  work. 


Short  Term  Releases 

During  the  Mountain  Iron  testing  series,  a  few  tests  were  con¬ 
ducted  with  release  times  considerably  shorter  than  the  standard. 
Although  such  tests  were  not  required  to  provide  the  final  equation 
sought,  additional  confidence  in  the  use  of  the  equation  seemed  possible  if 
such  data  were  available.  In  all,  9  five  minute  tests  were  made 
from  Source  A,  with  two  being  night  tests.  The  exposure  data  from  the 
5-minute  tests  formed  a  very  homogeneous  set,  with  relatively  little 
scatter.  Comparison  with  data  from  30-minute  tests  showed  that  the 
5-minute  tests  resulted  in  the  same  level  of  exposure  at  moderate  dis¬ 
tances,  1  to  3  miles  (Figure  5).  The  striking  difference  lies  in  the 
dependence  on  X  of  the  5-minute  tests:  whereas,  the  30-minute  tests 
decrease  with  the  -1.8  power  of  distance,  the  5-minute  tests  decrease 
with  the  -2.2  power  of  distance.  This  is  in  accordance  with  theory,  at 
least  in  principle,  if  not  in  magnitude.  The  30-minute  tests  are  little 
affected  by  longitudinal  diffusion,  but  the  5-minute  tests  must  be  sig¬ 
nificantly  affected,  since  the  longitudinal  dimension  of  the  short  tests  is 
usually  comparable  with  the  horizontal  dimension. 


Trajectory  Studies 


Several  techniques  for  estimating  the  trajectory  of  the  material 
released  during  Mountain  Iron  tests  were  investigated  to  determine  as 
accurately  as  possible  the  relative  usefulness  of  each. 

•  The  simplest  approach  involved  the  use  of  wind  direction  from 
a  single  WIND  sensor  to  predict  the  trajectory. 

•  A  more  complex  method  was  the  use  of  100-foot  winds  wherever 
available  and  54-foot-or  surface  winds  at  other  sites. 

•  Another  method  entailed  the  use  of  5-minute  averages  of  surface 
winds  at  all  possible  sites  to  construct  a  rather  detailed 
trajectory  prediction. 

These  methods  will  be  discussed  individually  to  provide  a  basis  for 
operational  trajectory  predictions. 
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Trajectories  from  30-Minute  Average  Winds 

The  possibility  of  using  rather  long-term  averages  of  wind  direction 
and  speed  to  estimate  plume  travel  was  investigated  by  using  30-minute 
average  data  from  the  operational  WIND  sites.  Only  a  few  tests  were  con¬ 
ducted  during  the  brief  interval  when  the  entire  WIND  system  was  up, 
prior  to  the  complete  shutdown  involved  in  the  installation  of  the  new  com¬ 
puter  during  June  and  July.  The  tests  were  MI-67,  and  69  through  96,  of 
which  14  tests  were  from  Source  A  and  5  from  Source  B. 

The  simplest  approach  to  trajectory  forecast  is  to  choose  a  single 
wind  sensor  to  indicate  the  wind  field  carrying  the  plume.  To  determine 
the  success  of  various  choices  of  wind  sensors,  the  difference  in  direction 
between  the  forecast  plume  direction  and  the  actual  directipn  was  calcu¬ 
lated  (in  degrees  azimuth)  for  all  the  tests  listed  above,  along  with  the 
standard  deviation  of  error  (Table  I).  Since  the  one -sigma  limit  is 
tabled  for  each  site,  the  trajectory  forecast  from  a  particular  site  fell 
within  the  listed  limits  of  error  68%  of  the  time. 

TABLE  1.  Plume  Centerline  Location  Error  Using  Selected 
Single  Wind  Observations 


Source 


Height  of 

Wind  Observation, 
_ ft _ 

Surface 

300  -  100 


Degrees  Error 
un  68%  of  cases) 


Surface 


Surface 

Surface 

Surface 


From  the  Table,  if  only  a  single  site  is  to  be  used,  it  is  clear  that 
Site  100  provides  the  best  prediction  of  plume  travel.  Thus,  for  a 
release  at  the  VIP-1  Pad,  a  centerline  drawn  according  to  Site  100 
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wind  direction  would  be  less  than  G  degrees  in  error  68%  of  the  time. 

A  similar  degree  of  success  is  plausible  for  releases  from  nearby  points 
(such  as  Pale  II),  especially  since  a  quarter  of  the  data  involv  'd  in  deriv¬ 
ing  Table  I  were  from  Source  B  releases,  rather  distant  for  the  VIP-1  pad. 

Nevertheless,  the  use  of  a  single  point  for  estimating  trajectories 
is  to  be  discouraged,  especially  in  such  terrain  as  found  at  South 
Vundotiberg,  since  the  probability  of  strong  deviations  in  direction  with 
increasing  distance  is  high.  The  80  minute  average  wind  conditions  from 
all  the  WIND  system  sensors,  using  100-foot  winds  where  available,  pro¬ 
vided  a  much  more  detailed  basis  on  which  to  work.  To  construct  the 
SO- minute  trajectories,  the  wind  field,  indicated  by  the  30-minute  average 
winds,  was  drawn,  and  the  path  of  a  parcel  starting  from  the  source  point 
was  drawn  according  to  the  field.  The  next  30  minute  average  field  was 
drawn,  the  path  of  the  parcel  for  the  second  30  minute  interval  deduced, 
and  drawn  onto  the  first  30  minute  map.  This  path  was  then  comoared 
with  the  observed  plume. 

Two  variations  were  investigated  at  this  time.  The  100  foot  wind 
field  was  drawn  both  with  and  without  the  information  received  from  the 
source  point  wind  sensor,  which  is  not  a  WIND  site  for  the  Source  A 
releases.  Likewise,  the  wind  fields  as  indicated  by  the  surface  winds 
were  drawn  with  and  without  the  source  point  wind.  The  results  are 
shown  in  Table  II;  although  the  lack  of  a  source  point  wind  is  a  serious 
handicap,  it  is  interesting  to  note  that  the  relative  success  of  this 
detailed  procedure  is  no  better  than  that  of  t.he  single-wind  process. 

For  example,  at  2  miles,  the  best  estimate  is  from  the  100-foot  winds, 
but  the  one -sigma  error  is  still  6  degrees. 

The  surface  winds  have  a  larger  error,  8  degrees,  which  seems 
to  Indicate  that  100  feet  is  a  better  estimate  of  the  transport  height  of  the 
bulk  of  the  plume.  There  are  no  definitive  data  on  this  point,  so  it  must 
be  conjectural.  It  was  shown  earlier  (in  the  section  on  aircraft  sampling) 
that  the  plume  regularly  was  sampled  at  hundreds  of  feet  above  the 
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TABLE  II.  Trajectory  Errors  in  Tests  MI- 07  and  MI- 09  to  *00 
Using  Selected  Wind  Data 

Pearces  Error  0H%  of  Cases 


Wind  Data 

1.8  mi 

3 .6  mi 

100  ft  height  with  source  wind 

JG 

■14 

100  ft  height  without  source  wind 

±12 

±10 

Surface  winds 

18 

I  1  3 

Surface  winds  without  source  wind 

+21 

±2  1 

terrain;  therefore,  a  great  deal  of  the  plume  was  well  removed  from  the 
surface,  at  least  during  daytime  conditions,  to  which  aircraft  sampling 
was  restricted. 

Trajectories  from  5-Minute  Average  Winds 

A  detailed  investigation  of  trajectories  as  indicated  by  the  surface 
winds  averaged  over  5-minute  intervals  was  also  carried  out  for  one  ease. 

The  wind  field  was  drawn  in  detail  from  all  available  wind  information  for 
each  five  minute  interval;  from  this  the  plume  travel  was  deduced  and  drawn  on 
a  map  and,  as  would  be  expected,  this  technique  yielded  the  best  estimate 
of  plume  travel.  The  results  of  the  trajectory  predicted  and  observed  for 
Mountain  Iron  Test  72  is  shown  in  Figure  G.  It  is  worthwhile  to  note  that 
the  curvature  and  spread  of  the  ground  pattern  of  exposure  is,  indeed,  indi¬ 
cated  in  the  wind  pattern  as  well.  No  attempt  was  made  to  relate  the  spread 
of  the  estimated  plume  segment  trajectories  to  the  observed  plume  spread, 
but  such  a  relationship  would  be  worthy  of  future  study. 

The  success  of  any  trajectory  analysis  probably  lues  in  accurate 
and  rapid  wind  information.  The  WIND  system  should  provide  an  adequate 
basis  for  well  estimated  plume  paths,  particularly  if  the  trajectories  are 
constructed  by  one  familiar  with  the  peculiarities  of  orographieally  dis¬ 
torted  wind  fields. 

One  of  the  by-products  of  the  trajectory  investigation  was  noting 
the  persistent  lack  of  correlation  between  the  average  wind  field  and  the 
flow  indicated  by  three  of  the  WIND  system  sensors:  010  in  Deal'  Creek 
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Canyon,  Oil  in  Spring  Canyon,  and  012  in  Honda  Canyon.  These  sites 
appear  to  be  very  strongly  influenced  by  the  nearby  walls  of  their  respective 
canyons,  so  that  their  information  is  not  at  all  indicative  of  the  flow  only 
a  few  hundreds  of  yards  away.  These  sensors  should  be  left  where  they  are. 
only  if  the  data  from  these  canyons  are  of  primary  importance.  If.'  their 
indications  are  used  uncritically  to  construct  wind  fields,  severe  and 
erroneous  curvature  in  the  wind  patterns  will  often  result. 
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CHAPTER  3. 

PHYSICAL  AND  CLIMATOLOGICAL  DESCRIPTION 
OF  SOUTH  VANDENBERG 

INTRODUCTION 

A  very  brief  description  of  the  geography  and  climatology  will  be 
given  here.  In  other  chapters,  the  relation  between  terrain  and  diffusion 
will  be  discussed.  The  relatively  short  time  available  for  collection  of 
imatological  data  precluded  any  thorough  climatological  summary. 
However,  the  climatology  published  by  Meteorology  Research,  Inc.  ,  in 
conjunction  with  an  earlier  study,  is  still  appropriate,  and  in  many 
places  provides  insight  into  the  meteorological  phenomena  at  South 
Vandenberg. 

THE  SITE 

South  Vandenberg  lies  on  a  westward  jutting  portion  of  the 
California  coast  about  a  hundred  miles  west-northwest  of  Los  Angeles 
(Figure  7).  The  Coast  Range  Mountains  dwindle  into  the  ocean  along  this 
section  of  the  coast,  and  form  the  eastward  buttresses  of  the  Santa  Ynez 
Mountains  which  culminate  in  4000-foot  high  ridges  several  miles  east 
of  South  Vandenberg.  The  mountains  are  characterized  by  ridges  1000 
to  1500  feet  high,  running  generally  east-west,  with  2000-foot  high 
peaks  along  the  southern  portion  of  South  Vandenberg.  The  most 
striking  feature  is  probably  Honda  Canyon,  a  rather  deep  and  steep¬ 
sided  canyon,  walled  on  the  south  by  Honda  Ridge  and  on  the  north  by 
Target  Ridge.  The  canyon  is  some  700  or  800  feet  deep  along  most  of 
its  length,  and  is  generally  no  more  than  1  1/2  miles  across.  Smaller 
canyons,  which  are  not  so  deep  but  even  narrower  (opening  mostly  north 
and  east),  branch  from  these  two  main  ridges.  The  result  is  a  broken 
appearance  which  is  much  more  rugged  than  the  description  in  feet  and 
miles  implies.  The  prominent  terrain  features  on  South  Vandenberg  are 
identified  in  Figure  8.  The  northwest  portion  of  South  Vandenberg  is  the 
Lompoc  Terrace --several  square  miles  of  level  and  rolling  grassland 
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FIGURE  8.  Prominent  Terrain  Features,  Mi-Phase  I 
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only  a  IVw  hundred  feet  above  the  ocean.  North  of  this  plateau  lies  the 
very  flat  flood  plain  of  the  old  Santa  Ynez  River;  from  this  smooth  river 
bottom,  the  Santa  Ynez  Mountains  rise  abruptly  like  a  pile  of  sand  on  a 
table  top.  The  coastline  lies  nearly  north- south  along  the  portion  of  South 
Vandenberg  north  of  Honda  Ridge,  then  turns  more  than  90  degrees  to 
east-west  along  the  southern  border  south  of  Point  Arguello.  This  abrupt 
curvature  of  the  coast  appears  to  be  related  to  the  common  occurrence  of 
northwest  winds  over  South  Vandenberg  that  veer  to  west  winds  over  Point 
Conception  a  few  miles  south  of  Point  Arguello, 

The  major  ridges  on  South  Vandenberg  are  mostly  covered  with  a 
dense  stand  of  low  growing  chapparal  and  related  shrubs,  as  shown  in 
Figure  9,  whereas  the  low  lands  are  generally  grassland- -the  floor  of 
Honda  Canyon,  the  Sudden  Ranch  area,  and  the  Lompoc  Terrace  being 
examples.  Trees  occur  in  scattered  clumps  only,  except  east  of  the 
head  of  Honda  Canyon  and  along  creek  beds.  The  vegetation  types  change 
rather  rapidly  with  increasing  distance  from  the  ocean;  trees  occur  in 
increasing  number  following  the  generally  increasing  average  precipita¬ 
tion  inland  in  this  area. 

CLIMATOLOGY 

The  period  of  time  during  which  Operation  Mountain  Iron  Phage  I 
was  fully  operational  spanned  only  eight  months,  much  too  short  for  a 
climatological  summary.  However,  wind  data  from  several  sites  were 
available  for  many  days  of  the  months  between  January  1966  and  July 
1966.  The  midseason  months-- January,  April  and  July--were  chosen 
for  complete  reduction  of  all  wind  data  available  from  several  sites: 

100,  VHF,  200,  Target  Ridge,  101,  055  and  301  (when  available).  The 
wind  data  were  read  as  30-minute  averages  centered  on  the  quarter 
hours,  with  direction  read  to  8  points  of  the  compass  and  speed  in  5 
knot  increments  or  calm. 

The  data  were  then  summarized  by  averaging  over  all  days  of  the 
month  and  over  a  given  time  interval.  The  time  intervals  chosen  were 


FIGURE  9.  Vegetation  on  Major  Ridges  of  South  Vandenberg 
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1/2,  1,  2,  3,  4,  and  6  hours,  starting  at  0000  hours.  A  tabulation  of  the 
summary  is  included  as  Appendix  B  in  Volume  I.  Included  in  the  monthly 
summary  for  each  site  is  the  highest  observed  30  minute  average  wind 
speed,  the  direction  from  which  it  came,  the  percentage  distribution  of 
wind  directions  in  the  given  time  period,  and  the  average  speed  during 
that  time  period. 

A  fairly  complete  picture  of  wind  patterns  over  Southern 

(4) 

California  is  given  by  Demarrais,  Holzworth  and  Hosier  for  the  mid¬ 
season  months.  By  use  of  their  regional  analysis  as  a  backdrop  and  the 
data  presented  here  for  details  of  flow  over  South  Vandenberg,  a  fairly 
reliable  generalized  picture  of  transport  near  Vandenberg  should  be 
possible.  To  illustrate  the  patterns  over  South  Vandenberg,  a  series  of 
maps  are  included  in  Appendix  C,  Volume  I;  these  maps  show  the  most 
probable  1-hour  average  wind  direction  for  day  and  night  conditions 
during  the  3  midseason  months. 
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MOUNTAIN  IRON  PHASE  I  TEST  PLUME  PATTERNS 


Some  of  the  113  tests  from  two  generation  (source)  points  on 
South  Vandenberg  are  sufficiently  interesting  to  include  for  reference. 
It  must  be  understood  that  the  plume  patterns  presented  are  not  defini¬ 
tive  lor  diffusion  prediction,  but  they  do  provide  an  orientation  and 
insight  into  the  cliniotological  regime. 
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APPENDIX  B 


MONTHLY  AVERAGE  SURFACE  WINDS 
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